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1. INTRODUCTION 



In relativistic heavy-ion collisions one is hoping to produce conditions where energy 
densities are high enough so that a new state of matter- the quark-gluon plasma can be 
produced. This state of matter lasts for a short period of time following the collision and 
may or may not be in equilibrium. Following this phase a transition to ordinary hadronic 
matter takes place and many of the processes which occur during the quark-gluon plasma 
phase might be masked by processes which occur in the hadronic phase. In order to 
determine processes which might be signals of the quark gluon plasma one needs to know 
the dynamical evolution of the plasma. This is because the particles that get produced 
during that phase have to travel through a time evolving plasma. In order to study this 
problem one needs a different way of thinking about field theory. Traditionally experiments 
in Elementary Particle Physics are black box experiments where initial particles enter a 
region, final particles exit the experimental region and all that is asked is how many 
particles of what type, energy, etc. enter various detectors. This type of experiment 
requires only a covariant S-Matrix theory to predict the probabilities to be expected in the 
detectors. However, if we want to know signatures of the quark gluon plasma,we actually 
need to follow the time evolution of the plasma and fields produced following the heavy 
ion collision. This requires a non-covariant real time formalism for the time evolution 
of the quantum fields. In these talks we would first like to discuss various formalisms 
for doing real time calculations in quantum field theory and then study in detail a very 
simplified model of the production of the quark-gluon plasma- Schwinger's mechanism for 
Pair-production from strong "classical" gauge fields. The value of doing a "first principles" 
calculation at this time, even if it is over-simplified, is multifold: 

(1) We can test the validity of existing semiclassical transport models of lepton production 
from the quark-gluon plasma. We have already discovered that these models have to 
be modified to correctly include Pauli blocking and Bose enhancement effects which 
were ignored. 

(2) We can determine the effective hydrodynamics and show that certain kinematic as- 
sumptions automatically lead to fiat rapidity distributions independent of the form of 
the equation of state. 

(3) We can determine the dynamical equation of state and in the next order in a systematic 
calculation in powers of (1/N) we will be able to study whether equilibration will occur 
and calculate self consistently lepton pair production rates. 
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First I would like to list the various approaches available to study real time processes 
in Quantum Field Theory. Each of these approaches needs an approximation scheme 
to reduce the number of degrees of freedom in order to make the problem numerically 
tractable. Three methods that my collaborators and I have studied in detail are: 

1- Functional Schrodinger Equation + variational approximations fll] 

2- Truncated Heisenberg Equations - Large-N expansion or Mean-field approximations 
to the Dyson equations [g 00 0. 

3- Schwinger's closed time path - Path Integral Formalism in a large-N expansion 

There also exists an alternative formalism related to the truncated Heisenberg equa- 
tions based on the Wigner Distribution function which has been discussed by Rafelski and 
his collaborators [|T^. Once we have chosen a method we have to decide how to specify 
the initial data at t=0. In these different approaches we have to specify 

1- Initial position and width of say a Gaussian Wave Packet at t=0 in the Schrodinger 
picture 

2- Number and pair densities etc. in Heisenberg picture. 

3- Initial density matrix in the Path Integral Approach. 

In classical field theory, such as classical electrodynamics, the theory is finite and any 
smooth initial configuration of the field is allowed for the initial value problem. When we 
have a semiclassical field theory for the expectation value of the fields however, the initial 
data can be reinterpreted in terms of the particle language and even a smooth initial 
configuration of the field might not be consistent with certain physical constraints such as 
the initial state having finite number density at t=0 with respect to an adiabatic vacuum. 
(This requirement is automatic for finite temperature field theory). Thus arbitrary initial 
data may not be consistent with renormalizability. This is discussed in detail in 
We also have an additional new problem to face - how to perform renormalization in a 
non-covariant formulation of the field theory. To do this we isolate the divergences in an 
adiabatic (WKB) expansion of Green's functions. This method is similar to the technique 



of adiabatic regularization used by Parker and Fulling [|T3[ in their study of semiclassical 
gravity. The problem we will address in detail in these lectures is pair production of either 
Bosons or Fermions from strong Classical Fields which are either functions of time t, or 
fluid proper time r = (t^ — 2^)^/^. We will compare the results of the numerical simulation 
of this problem (for the degradation of the field, the particle spectra, etc.) with a semi- 
classical transport approach using a Schwinger-inspired source term [|14[[0|jT^. We will 
also discuss the effective hydrodynamics derived from the expectation value of the energy 
momentum tensor of the quantum theory. 
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2. SUMMARY OF THE DIFFERENT STRATEGIES IN X^p^ FIELD THE- 
ORY 

For simplicity let us first study these different approaches to initial value problems in 
the simplest case- Xip^ field theory, 
a) Schrodinger Picture: In the Schrodinger picture the Initial State is described by a 
wave functional at t=0. For example a Gaussian wave functional is 

= exp[-/" [ipix)-0ix)][G-\x,y)/4~iJ:ix,y)]My)-0{y)]] (2.1) 

Jx,y 

The time evolution is given by the Functional Schrodinger equation |jl[]: 

idip/dt = Hip 

H = J d'x[~5y6^' + ^{Wcpf + V{cp)] (2.2) 

This is a generalization of the usual Schrodinger equation: 

tp{x) =< tp\x >, X ^ ip{x, t); 
p = —id/dx ^ TT = —id/d(f 

idt(;/dt = Ht(;;H = -d^/dx"^ + V{x) (2.3) 

with initial condition: 

?/^(0) = exp[-a(x -Xo)2]. (2.4) 

One might imagine solving (2.2) on a computer by introducing a lattice in d dimensions 
and converting the functional derivatives into partial derivatives. One then quickly realizes 
that the number of degrees of freedom in equation (2.2) is rather overwhelming. To control 
this problem one uses variational trial wave functionals which become "exact" in the large- 
N limit- namely Gaussians. The equations of motion for the variational parameters can 
be obtained from Dirac's variational principle ||17|| : 



T = j dt<^\id/dt-H\^ > (2.5) 
5r = ^ Schrodinger's equation: 

id/dt- H\-<l! >=Q (2.6) 
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In the representation one can choose a Gaussian trial wave functional: 

< <^|*„ >= ipv['P,t] = exp[- / [ip{x) - (p{x,t)] 

Jx,y (2.7) 

[G~^{x, y, t)/4 - iE{x, y, t)] [ifi{y) - ifi{y, t)]] + m{x, t) [ip{x) - ifi{x, t)] 
where the variational parameters have the meaning: 

G{x,y,t) =< ^y\ip{x)ip{y)\^y > -ip{x,t)ip{y,t) (2.8) 
Then the effective action for the trial wave functional is 
T{ip, TT, G, S) = J dt< ^y\id/dt - H\-^y > 

= J dtdx[7r{x,t)d(p{x,t)/dt+ J dtdxdyT,{x,y)dG{x,y,t)/dt 

-J dt<H> (2.9) 



where 



< 



H>= j dx{7rV2+2EG'E+G'-V8 + V2(V^)^-l/2V^G'+l/2F"[^]G'+l/8F""[^]G^}. 



< > is a constant of the motion and is a first integral of the motion. For \(p^ field 
theory we get the following equations of motion: 

TT{x,t) = V'^ip-d <V > /dip; 

(fi{x, t) = TT 

G{x,t)^2 J dz[J:{x,z)G{z,x) +G{x,z)J:{z,x)] 

t{x,t) = -2 J dz[^{x,z)^{z,x) + G-'^/8 

+ [l^l-d<V> /dG]S%x-y) (2.10) 

If there is translational invariance and ip=0 we obtain a second order differential 
equation for G(k,t), the Fourier transform of G(x,t): 

2G{k, t)G{k, t) - G^{k, t) + AT{k, t)G^{k, t) - 1 = 

T{k, t)=k^+ m'^{t);m^{t) = fi^ + ^X J [dk]G{k, t) (2.11) 



This approximation is called the time-dependent Hartree-Fock Approximation and is 
equivalent to the leading term in a 1/N expansion of the field theory 0. To understand 
this trial wave function let us look at a simple quantum mechanics problem- the harmonic 
oscillator with a gaussian initial state. Harmonic oscillator: V(x) =1/2 m x^, 

Initial conditions: 

^'(x,0) = [27rG'(0)]-^/2exp{-a;V[4G'(0)]} 

q{0) =<x>=0 (2.12) 

For the harmonic oscillator a Gaussian remains Gaussian as time evolves so that 

= {2nG{t))-^/^exp{-x^[G-\t)/4 - zE(t)]} (2.13) 

We find that the conserved Energy can be written in terms of G as follows: 

E=<H>= G^/8G + Gm^/2 + G'^/S = G'^/SG + V[g] (2.14) 

We plot V[g] in fig 1. From fig. 1 we see that the ground state is G = l/(2m). If at 
t=0, Go = 1/(2M) ; m M then 

G{t) = 1/2(^0 + Gi) + 1/2(^0 - Gi) cos(2m(t - to)) (2.15) 

Thus the width oscillates with frequency 2m between Go and Gi. Generalizing to 
free field theory (which is just independent harmonic oscillators) we have instead for each 
mode of momentum k: 

< H{k) >= G^/8G + {k^ + m^)G/2 + G'^/S (2.16) 

This leads to the same result for G(k,t) as for G(t) with m— ujk = (/c^-fm^)^/^. However in 
field theory, unlike quantum mechanics, an arbitrary initial Gaussian state is not necessarily 
a physically valid choice since it might correspond to an infinite particle density or energy 
density when compared to the adiabatic vacuum. Thus the particle interpretation implies 
that one needs to restrict the large k behavior of G(k) at t=0 to be a physically allowed 
initial state with finite particle number, energy density etc. Otherwise one gets extra 
unwanted infinities in loops, 
b) Heisenberg Picture: Green's function approach 
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In problems where there is spatial homogeneity one has a Fourier decomposition for 
a charged field <^ in terms of mode functions fk (t) which depend only on the time and the 
usual creation and annihilation operators a and b which satisfy the canonical commutation 
relations: 

^x,t) = I [dk][fk{t)ake"^^ + f:{t)hte-'^-] 

K,a+] = = {27:)H\k-k') (2.17) 

The initial state |i > is totally specified by specifying at t=0 the matrix elements of a and 
b: 

< i\alak\i > = {2TiY8'^{k - k')n+{k) 

< i\bkak\i > = {2T^Y5'^{k + k')F{k)etc. (2.18) 
The equation for the expectation value of the equation of motion is: 

<i\{-B + m^)^ + X{^+^)^\i>^0 (2.19) 

We see from these equations that we also need to solve the equation of motion for 

< i\X{ip^ip)ip\i >. 

In general we get a Heirarchy of Green's function equations- The BBGKY heirarchy. 
To make practical progress we need a truncation scheme which allows us to solve the lowest 
order problem and then systematically calculate corrections. In the large N expansion the 
lowest order approximation leads to a factorization 

< i\{ip~^(p)(p\i > —< i\{ip^ip)\i >< i\(p\i > 

= G{x,x;t) <i\(p\i> (2.20) 

where the fourier transform G(k,t) of G(x-x'; t) obeys the same equation as the width of 
the Gaussian wave function in the Schrodinger equation in the Hartree approximation. 

2G{k, t)G{k, t) - G^{k, t) + 4r{k, t)G^{k, t)-l = 
T{k,t) = k'^ + m^{t);m'^{t) = fj,^ + J [dk]G{k,t) 

G{x,x;t) = J[dk]G{k,t) (2.21) 



Thus the large-N expansion (Hartree approximation, mean field approximation) truncates 
the hierarchy of coupled Green's function equations making it necessary to only solve the 
coupled one and two-point Green's function equations. 

In these mean field equations the problem reduces to an external field problem in that 
the quantum field <^ obeys the equation: 

{-B + rr?{t))<fi = (2.22) 

Because we have an external field problem with spatial homogeneity: the mode functions 
f(t) in (2.17) obey: 

(a^ + a;^)/ = 0; a;^ = /c^ + m^(t) (2.23) 
The canonical commutation relations lead to a constraint on the mode functions: 

fktk - fkfk = ^ (2.24) 
which is automatically satified by the WKB form ansatz: 

/fc(t) = [20K(t)]-'/'expH2/fe(t)] 

ykit)-^k{t) (2.25) 

which lead to the equation 

nlit) + nk/{2nk) - ^{nk/nkf = u;l{t). (2.26) 

At t=0 one has in general for the initial state: 

< i\a^ak\i >= {2TTf5'^{k - k')n+{k) 

< i\bkak\i >= (27r)^5^(A; - k')F{k) 

For an adiabatic vacuum: n(k) =F(k) = 0, and the initial conditions on 0, are 

i}(k, t = 0)= uj{k, t = 0); (l{k, t = 0) = Cj{k, t = 0). (2.27) 

This formalism, however is perfectly general and one could take any initial state with 
an integrable phase space particle density n(k) and pair density F(k). As a particular 
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choice one could have chosen at t=0 an equihbrium configuration of pions described by a 
temperature kT = (3~^ 

n{k) = l/{exp[PE{k)] - 1) (2.28) 

c) Path Integral Approach: Closed time-path formalism 

The only formalism that allows a systematic approach to initial value problems is the 
closed time-path approach of J. Schwinger|P] which was further elaborated by Keldysh^and 



put into a Path Integral framework by Chou, Su, Hao and Yu[^0[. This Path Integral ap- 
proach allows standard Path Integral approximation schemes such as the large N approx- 
imation as well as ensuring causality for the Green's functions for initial value problems 
||18|| . The starting point for determining the Green's functions of the initial value problem 
is the generating Functional: 



Z[J+,J ,p] =< z|T*(exp{- J iJ-ip-})\out X out\T{ex.p J iJ+ip+)\i > (2. 



29) 



This can be written as the product of an ordinary Path integral times a complex 
conjugate one or as a matrix Path integral. 

Z[J+, J~, p] = J dip~^dip~ < i\p\(p-, i > expi[{S[ip+] + J+(p+) - {S*[ip-] + J-^P-)] 

= j dipc,exv^{SYPa\ + Ja'^a) < ^i,ApW2,i> (2.30) 

where < ipj^i\p\(p-^i > is the density matrix defining the initial state. 
This leads to the following matrix Green's functions [|TT|: 

G++ - 6^ In Z/6J+6J+\j=o =< T{ip{xi),ip{x2) > 

In Z/6J-6J-\j=o =< T*{ip{xi),ip{x2) > 
G+- = In Z/dJ+dJ~\j^o =< (p{x2), i^ixi) > 

G-+ = 5^ In Z/6J-6J+\j^o =< ip{xi),ip{x2) > (2.31) 

The matrix Green's function structure insures causality. In this approach it is easy 
to generate a 1/N expansion in analogy with ordinary field theory. The diagrams are the 
same as in the usual 1/N expansion, except the Green's functions are the matrix Green's 
functions described above. If in lowest order in (1/N) we have an external field problem 
as described above, one can directly use the mode solutions of the previous methods to 
determine the lowest order matrix Green's function of eq. (2.31). This obviates the need 
to discuss the initial density matrix of the theory, since it is these Green's functions which 
then enter the diagrams of the higher order calculations. 
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3. MAIN EXPANSION IDEA: FLAVOR SU(N) 

In many problems one of the fields can be treated classically to first approximation- 
pair production in Strong Electric or Gravitational fields. This makes the lowest order 
problem an external field problem. One way to generate a systematic expansion whose 
lowest order is an external field problem is by introducing N copies of the original problem 
and expanding in Flavor SU(N). This is most easily done in the Path Integral formalism. 
For the initial value problem one would use the matrix Green's functions discussed above. 
Having an extra large parameter N allows an evaluation of the Path integral by Laplace's 
method (or the method of Steepest Descent). To obtain the large N expansion one realizes 
that if there are N flavors the loops carry an extra N. Rescaling the fields then display an 
overall factor of N in the effective action which includes the loops. Examples: 

A(^^ : X = ^'^ 

Z = J dxj d(fexp[- j {d^ipf + Axv'^ - Ax' + /iV + Jv + Sx] 
^ (^i, z = 1, 2, • • • AT; A ^ \/N- ^ iV^/Vi; X ^ ^^X, Ax ^ Ax- 

Integrating over we obtain: 

Z = Jdxexp{-N[x^ + ^TrlnG-^-jGj]} 

^ j c^Xexp{-iVS'e//(x)} 

G-^ = [-U + ^i^ + Xx]5{x-y) (3.1) 

Evaluating the Path Integral at the Saddle point, SSeff{x)/^X — leads to the self 
consistent external field problem 

[-□ + fj,^ + Xx]v = 0;x = v'' + G{xx) (3.3) 

In QED we obtain an external field problem by integrating out the fermions (which 
have now N flavors to give an extra N to the determinant) and then rescaling the flelds to 
display the overall factor of N in the effective Action: QED: 

Z = y dAf^ J d^d^ exp[ j dxi—F"^ + ^{i-^d - e /I + m)^'} + Wry + ry*] 
*^*i;e^e/yiV,^^^ViV (3.3) 



Integrate out the N species of fermions 
j dA^exp{-NS,ff{A^)} 

Seff{A^) = j dx^F^ + Trln{S-\x,y;A))+fjS{x,y;A)rj] 

S-\x,y;A) = {i'yd - e /A{x) + m)6{x - y) (3.4) 

Evaluating the Path Integral at the saddle point, dSeff{A^)/dA^ = leads to the 
external field problem: 

(i'yd - e /A + m)-^ = (3.6) 

where A is an external field, \E' is a quantum field. We also obtain the semiclassical Maxwell 
Equation: 

d^Ff"^ =< / >= e < ■*"7^^' > . (3.7) 

In all these problems one has in leading order in 1/N a straightforward problem of a 
quantum field theory in a background field which allows a normal mode decomposition in 
terms of the solutions of the classical field equations. Renormalization can be carried out 
by an adiabatic expansion of the mode equationjjl^. The effect of quantum fiuctuations 
about the semiclassical field can be systematically taken into account by calculating the 
fiuctuations about the leading stationary phase point in the Path Integral order by order 
in the 1/N expansion. 



4. PARTICLE PRODUCTION IN THE CENTRAL RAPIDITY REGION IN 
HEAVY ION COLLISIONS 

A popular picture of high-energy heavy ion collisions begins with the creation of a fiixx 
tube containing a strong color electric field[0. The field energy is converted into particles 
as qq pairs and gluons which are created by tunnelling- the so-called Schwinger mechanism 



|j2^ [|21[ [g2[ . The particle production can be modeled as an inside-outside cascade which is 
symmetric under longitudinal boosts and thus produces a plateau in the particle rapidity 
distribution. The boost invariant dynamics, in a hydrodynamical picture gets translated 
into energy densities (such as E^ ) being functions of the proper time. We take this as 
an initial condition on the fields in an initial value problem based on this pair-production 
mechanism. First let us look at the case where the electric field is a function of real time 
t, treating later the more realistic case where E= E we first 
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consider particle production from a spatially uniform electric field such as that produced 
between two parallel plates. This is an idealized model of a flux tube for QCD. The problem 
of pair production from a constant Electric field (ignoring the back reaction) was studied 
by J. Schwinger in 1951 PUI. The physics is as follows: One imagines an electron bound 



by a potential well of order |Vo| ~ 2m and submitted to an additional electric potential 
eEx (as shown in fig. 2 ). The ionization probability is proportional to the WKB barrier 
penetration factor: 

exp[-2 / dx{2m{Vo-\eE\x)y/^] =exp{--m'^/\eE\) (4.1) 
Jo 3 

A direct calculation due to Schwinger from first principles using the effective action 

in an arbitrary constant electric field (ignoring the back reaction) gives instead 

w = [aEV(27r2)]E^^i^ ^ exp(-n7rmV|e£;|). (4.2) 

This equation tells us that pair production is exponentially suppressed unless eE 
> Trm^. So we expect (as we find in fig. 3) that there is a crossover value of E where 
the time it takes for E to first reach zero (remember there are plasma oscillations) is 
relatively short. Schwinger 's result only applies when we can ignore dynamical photons 
(as well as back reaction)and is related to the lowest order in 1/N calculation where the 
electric field is treated as a classical object. Schwinger's analytical result was subsequently 
used as source term for an approximate transport theory [|T3[, []T^ approach to the 
back reaction connected with pair production which we will later compare with our exact 
numerical results. 

We will choose the electric field in the z direction and choose a particularly simple 
gauge: 

E = E{t)k ~A = A{t)k; E{t) = -dA/dt (4.3) 
To maintain spatial homogeneity we have from Maxwell's equation: 

V-E = p (4.4) 

that the plasma of produced particles must be neutral. In scalar QED, the equation for 
the quantum field ip is 

-(a^ -zeAa)(a" -zeA")$ + ;u2$ = (4.5) 
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and for the electromagnetic field: 



a^^F^" =< C{-ie(^+d^^ - - 2e2yl^$+$} > (4.6) 

where C denotes charge symmetrization with respect to and For our constraints on 
the field E and our choice of gauge we get: 

-dE/dt =< jz >=^J [dk] {kz - eA{t))G{k, t) (4.7a) 
where G{k, t) = [< (y?V + > -2(p*(p]FT (4-76) 

For QED we have instead the field equation: 

[i-fd - e^{t) -m]'i/{x,t) = (4.8) 

and the semiclassical Maxwell equation: 

-dE/dt =< jz >= < i\[^{x, t), 73*(x, t)] \i > (4.9) 

The fact that the external field is independent of space (spatial homogeneity) means 
that one has a simple normal mode expansion of the fields just as in Xip'^ field theory 
described earlier. 

For Scalar QED we have 

^x,t) = J [dk][fk{t)ake"^- + f:{t)hle-'^-] 

ul{t) ^ [k - eA{t)f + I? + ki (4.10) 
Repeating the arguments of (2.22 - 2.25) we again obtain for the generalized frequency 

Qlit) + Qk/i^Qk) - li^k/^kf = ujlit). (4.11) 
where now u) is given by (3.10) Spatial homogeneity requires translational invariance, 

W{x - x', t, t') = J [dk]W{k, t, t')e'^^''-'''\ 
12 



This in turn requires that 



< a^ak > 



< bkttk > 



{2T:f5\k + k')F{k) 



(4.12) 



Thus we obtain for G(k,t) 



G{k]t) = VL~^{k,i){l+n+{k)+n-{k) + 2F{k) cos[2yk{t)\} 



(4.13) 



This is the most general form of the propagator that one would use in the diagrams of 
the 1/N expansion, where n and F are the particle and pair phase space densities at t=0. 
These parameters also totally specify (in leading order in 1/N) the density matrix at t=0. 
To solve the field theory in leading order in 1/N (ignoring questions of renormalization 
to be discussed below) one solves the second order differential equation for each mode 
function ^k{^)i determines G(k,t) and then solves the back reaction equation: 




(4.14) 



For QED one has to deal with the spinor structure: 




(4.15) 



If we choose a basis where 7°7'^ is diagonal: 



7Vx. = KXs.s^ 1,2 ^ A = l;s = 3,4^ A = -1 



(4.16) 



Then the spinors u and v obey the equation 




(4.17) 



Squaring the Dirac equation by letting: 




(4.18) 
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we find that the mode functions / now obey: 



[d^+ul{t)-iXs7r]fk{t) = 0, 

n^k-eA (4.19) 
If the operators ak and bk obey the usual anticommutation relations: 

{aks, 4, J = {bks, 4, J = {2nf6^k - k')5ss' (4.20) 
the fk are constrained to satisfy 

^2j*aj/3 ^ j*ajP ^ ■^^^^j^*aj(3 _ y*«j/3] = §a(3 ^2 (4.21) 

Parametrizing the positive and negative frequency solutions: 

/±(t) = iV±exp / g±{T)dT, (4.22) 

we find: 

g+ = -[Xsn + Q]/2n-in (4.23) 
where the generalized frequencies, flkit) now satisfies the equation: 

nl{t) + VLk/{2VLk) - h^tlu/VLkf - 7rV(4f]2) - A.TrJi/ll^ (4.24) 

Ignoring renormalization, the solution of QED is obtained by simutaneously solving 
for these modes and also for E(t) which is obtained from the Maxwell equation: 

dE/dt = 2eS^^i j [dk]{pl + m2)A,|/+ (t)|2 (4.25) 

5. RENORMALIZATION 

The equations of the previous section as they stand are not finite in the continuum 
since the sum over modes in (4.14) and (4.25) contains a divergence related to the renor- 
malization of the charge (as well as the wave function) resulting from the charged particle 
loops in the definition of the current. 

Let us first look at Scalar QED where the back-reaction equation is: 
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-dE/dt =< j >= e J[dk]{k^ - eA{t))n-^[l + N{k)...] (5.1) 

We first see that N(k) has to faU fast enough at large k to not lead to any further 
divergences- this is equivalent to the condition that the initial number density p is finite. 
The integral of contains a divergence proportional to dE/dt which renormalizes the 
charge (as well as the field E). To isolate this divergence one makes an adiabatic expansion 
of the equation for the generalized frequencies Qk- That is, we imagine that the time 
derivatives are small d/dt^ e d/dt : 

e^[Qk/i2nk) - ^{Qk/^kf] = ujlit) - Qlit) (5.2) 
and we then expand in powers of e 

l/Qk = + e^cvk/^ujk - ^{uk/ujkf} + 0{e^oj^^)] (5.3) 

We see that terms with higher derivatives are associated with more convergence factors 
of 1/k so that one only has to consider the first two terms in the adiabatic expansion to 
isolate the divergences which are interpreted as the standard charge renormalization. The 
log divergence comes from the term 

ivk = e{dE/dt){k - eA)oj~^ (5.4) 

After integrating over k this leads to a term of the form: 

e^Se^dE/dt; <^e^ = ^ y [dk]uj^^ = 7r(0) (5.5) 

where 7r(0) is the usual vacuum polarization at q^=0. Subtracting this term from both 
sides of eq. (5.1) we obtain: 

-edE/dt{l + 6^77(0)) = e'^[J [dk]{k, - eA{t))G - eTT{0)dE/dt]. (5.6) 
The Ward identity tells us that eE — crEji; and the renormalized charge is determined 

by 

e| = eV(l + eV(0)) (5.7) 
so the explicity mode by mode finite renormalized equation is 

-dEn/dt = eR j [dk]{k - eA{t))[n-^ - lu'^ - e|(A; - eA{t)){dE/dt)io-yA] (5.8) 
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For QED one gets instead after charge renormalization: 

dEn/dt = 2ei?E^^i j [dk][{pi + m')\s\ft? ~ e^dEn/dt uj-^ (5.9) 

6. HEAVY ION COLLISIONS AND BOOST INVARIANT DYNAMICS 

In e"*" e~ annihilation, hadronic coUisions and in heavy-ion colhsions particle produc- 
tion in the central rapidity region can be modeled as an inside-outside cascade which is 
symmetric under longitudinal boosts which leads to a plateau in the particle rapidity dis- 
tributions. This boost invariance also emerges dynamically in Landau's hydrodynamical 



model |]2^ and forms an essential kinematic ingredient in the analyses of Cooper, Frye and 



Schonberg as well as Bjorken|2^. It was recognized by Cooper and collaborators and 
further elaborated by Bjorken that in a hydrodynamical framework scale invariant initial 
conditions : 

V = z/t,e{x,t) ^ e{T),T^ = t^ - (6.1) 

would automatically lead to fiat rapidity distributions. In the context of transport or field 
theory modelling of the heavy ion collision, after an initial time tq, energy densities are 
expected to be functions only of the fluid "proper time" r. We therefore assume that 
the kinematics makes the electric field E only a function of the proper time r. For this 
kinematical choice it is convenient to introduce new variables r, ry the fluid "proper time" 
and the fluid rapidity (when v=z/t) via : 

z = T sinh rj^t = T cosh rj. (6-2) 

This change of coordinates to (r, 77) from (t,z) can be accommodated by the usual formalism 



of curved space (except the curvature here is zero) . One introduces the metric in 

curved space 

^7c./3 = diag(-l,0,0,r2). (6.3) 

Maxwell's equations 

{-9)-"''dp[{-g)-''''F'^^]=3^ (6.4) 
becomes for an electric field E{t) in the z direction 

E{r) = F,t = Fr,r/T = -T-'drA{r) (6.5) 

-l/TdAl/rdrA{r)]=<r > (6.6) 
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For Scalar Electrodynamics the equation for x = V^V is 

{dl + r-\d, - ieA{r)f + 1/4] - dl - + m^)x = (6.7) 

The rescaled field x has the same Fourier decomposition as (f) had in flat space with 
the mode functions f obeying 

[d'r + u^l{r)]Mr) = (6.8) 



however now 



r) = [k- eA{T)YlT' + fci + + l/(4r^) (6.9) 



so that the longitudinal momenta get suppressed at large r. For fermions one has the added 
complication that the covariant derivative now has a spin piece: (denote the Minkowski 
indices with ct, /3 the curvilinear coordinates with \xv) 

= 9^ + - ieA^, 

E-^ = l[7",7^] (6.10) 
and the vierbein represents the transformation to the Minkowski coordinates: 

9y^. = V^V^^Vcp: t = I'^V^ (6.11) 

Maxwell's equation becomes: 

-T~^dE{T)/dT =< f > = -^e < i|[^',7 > 

= < ^|[^'^7V*]K > (6.12) 
zr 



and the fermion mode functions now obey 



[d^,+ul{T)-iX,ir]fk{T)=0 (6.13) 

where 

loUt) = n^+pl + ^^;n=ip- eA{T))/T (6.14) 

The divergences in Maxwell's equation in curved space can be renormalized as before 
by an adiabatic expansion in the variable r. The details of this calculation are presented 



m 
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7. PARTICLE PRODUCTION RATES AND THE BOGOLIUBOV TRANS- 
FORMATION 

The wave functions of the first order adiabatic expansion g^^^ f^{t) where 

flit) = {2ujk)-'/'eM-i f^k{t')dt'] (7.1) 

Jo 

form an alternative basis for expanding the scalar fields and allows one to define an in- 
terpolating number density N(k,t) which becomes the true one as t— > oo. Expanding the 
field in terms of f^(t) we have 



t) = / [dk]K{t)e^''^fl{t) + r_l{t)bl\t)e-^'^-]where alit ^ oo) = etc. (7.2) 



In this expansion the creation and annihiliation operators are time dependent. We also 
have our previous expansion in terms of the time independent operators a and b related 
to the initial state: 

^{x,t) = J [dk][Ut)ake^'^ + fl,{t)ble-^''% (7.3) 

We recognize that and a^t) are related by a unitary transformation. The Bogoli- 
ubov coefficients are defined by 

al{t)=a{k,t)ak + f3*ik,t)bl 
bl{t) = a{k,t)bk + P*{k,t)al 
\aik,t)\' + \Pik,t)\' = l (7.4) 

The number of particles produced per unit volume is just 

V-^dN/dk =< t = 0|6°"* +6f ^ + a°"^ ^a^"*|t = > (7.5) 

The interpolating number density is defined in terms of the first order adiabatic op- 
erators: 

V-'dN{k,t)/d^k = {t = 0\bl\t)bl{t) + al\t)aUt)\t = 0) 

= (1 + N(k))\/3\'^ + N{k)\a\'^ + 2Re{aPF{k)} (7.6) 
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For N=F=0 (the adiabatic vacuum at t=0) 

v-^dN{k,t)/dk = = {4uknk)-^[{nk-uJkf + ^{nt/nk-dJk/ookf] (7.7) 

We see that adiabatic initial conditions (no particle production at t=0) are 



Ofc(O) =cc;fc;Ofc(0) = Wfc(O) (7.8) 
For fermions we have instead: 

V-^dN{k, t)/dk = T.M^ + A7r)(2u;)-i[wV+l' + - ^cc;(/*+c>o/+ - /+ao/*+)]. 
Similar expressions exist for the boost invariant problem . 



8. TRANSPORT APPROACH TO MULTIPARTICLE PRODUCTION 

A classical kinetic theory approach to the back-reaction problem as discussed in 



[1T5[| [|T6| introduces a phase space single particle distribution function f{x,p,t) in the 
presence of a homogeneous electric field and with a phenomenological source term inspired 
by Schwinger's solution for the constant external field. 

df/dt = df/dt + eE(t)df/dp 
= dN/dtdzdp 

= \eE{t)\\n[l ± exp[-nm'i/\eE{t)\]]6{p) (8.1) 

± stand for boson(+) or fermion case (-). The right hand side of (8.1) is a naive use 
of Schwinger's formula (valid when no particles are present and for constant fields with E 
replaced by E(t) {or E(r)}. This approach was recently used to predict dilepton production 
from the quark-gluon plasma |2^]. A potential problem with replacing constant E by E(t) 



is that in the field theory simulations E(t) is rapidly varying in time. A more serious 
problem is that once particles are produced, Schwinger's derivation, which was for particle 
production from the vacuum, is no longer valid. This however can be fixed up by the 
following arguement. Once particles are present there is an additional quantum mechanical 
effect due to statistics- Bose enhancement or Pauli Suppression. For the external field 
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problem one always has a normal mode decomposition at each time t. Thus the creation 
and annihilation operators at different t are again connected by a unitary transformation: 



Therefore 



b{k, t + At) = a{t + At)b{k, t) + [3{t + At)S{k, t) 
lap + = 1; \h^h\ = n+; \a^a\ = n-;n+ = ri- = n (8.2) 

n{t + At) = n{t) + 2\(3{t + At) pjl ± n} (8.2) 



or 



An/At = 2|/3|2{1 ±n}/At (8.3) 

where the +(-) is Bose enhancement (Pauli suppression). The Pauli suppression ensures 
n(k) < 1 for fermions. Thus to include this effect we will modify the right hand side of 
(8.1) by multiplying by (1± 2 f (p,t)). This modified transport eqaution, as we will show 
below gives much better agreement with the field theory calculation. One can solve the 
Viasov equation using the method of characteristics. From dp/dt = eE and f(p,0)=0 one 
obtains: 

/(p,t) = E,ln[l±exp[-7rmV|eE(tO|]] (8.4) 
where the t^ are determined from 

p + eA{t)+eA{ti) = 0]ti<t (8.5) 

The back reaction equation is now 

d^A/dt^ = jcond + ipol 

icond = 2e J[dp]pf{p,t)/{p^ + m'')^/^ 

jpoi = 2/E J [dp]{p'^ + m^f/'^d^N/dtdxdp (8.7) 

where 

dN/dtdxdp = 

(1 ± 2/(p, t))\eE{t) I ln[l ± exp[-7rml/\eE{t) \]]d{p) 
A similar expression holds in boost invariant dynamics as discussed in . The transport 



approach with the enhancement (suppression) factor gives reasonable agreement with the 
direct numerical solution of the field theory (in lowest order in 1 /N) as long as we coarse 
grain the field theory result in momentum bins. 
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9. HYDRODYNAMIC CONSIDERATIONS: ENERGY FLOW 

From a hydro dynamical point of view, fiat rapidity distributions seen in multiparticle 
production in p-p as well as A-p and A-A collisions are a result of the hydrodynamics 
being in a scaling regime for the longitudinal flow. 

That is for v=z/t (no size scale in the longitudinal dimension) the light cone variables 

2; = Tsinh?/; t = r cosh?/ (9-1) 
become the fluid proper time t = t(l — v^)-^/^ and fluid rapidity: 

T] = 1/2 ln[(t - x)/(t + x)] 1/2 ln[(l - v)/{l + v)] = a (9.2) 

In the rest frame (comoving frame) of a perfect relativistic fluid the stress tensor has 
the form: 

= diagonal (e, p, p, p) (9.3) 
Boosting by the relativistic fluid velocity four vector vJ^{x,t) one has: 



{e + p)u^u'' -pg^"" (9.4) 



Letting vP = cosh a; vP = sinhct, we have when v = zt that rj = a, the fluid rapidity. 
If one has an effective equation of state p = p{e) (which happens if both p and e are 
functions of the single variable (r) as well as for the case of local thermal equilibrium) 
then one can formally deflne temperature and pressure as follows: 



e + p = Ts;de = Tds;lns^ Jde/{e + p) (9.5) 

Then the equation: 

u^d'-T^^ = 

becomes: 

a^(s(T)«,) = (9.6) 

Which in 1 + 1 dimensions becomes 

(is/cZr + s/r = or ST = constant (9.7) 

The assumption of hydrodynamical models is that the initial energy density for the 
flow can be related to the center of mass energy and a given volume (say of a Lorentz 
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contracted disk of matter). It is also assumed that the flow of energy is unaffected by the 
hadronization process and that the fluid rapidity can be identifled in the out regime with 
particle rapidity. Thus after hadronization the number of pions found in a bin of fluid 
rapidity can be obtained from the energy in a bin of rapidity by dividing by the energy 
of a single pion having that rapidty. That is one assumes that when the comoving energy 
density become of the order of one pion/(compton wave length) |2^] ||2^ we are in the out 
regime. This determines a surface deflned by 

e(r/) = m^/K (9.8) 

On that surface of constant r 

dafj, = ATXo?{dz, —dt) = An t f {cosh rj, — sinhry) 
dN/drj = Ana? /mT^[{e -\- p) coshacosh(?7 — a)pcoshrj\/ cosha = Ano?'(i{Tf) /m-n: (9.9) 

which shows that when rj = a one gets a flat distribution in fluid rapidity. An extra 
assumption is needed to identify fluid rapidity a with particle rapidity y = l/21n[(i?7r + 
Pt^)/{Et^ — Ptj:)], where Pt^ is the longitudinal momentum of the pion. 

What I would like to show next is that in a fleld theory calculation based on the 
Schwinger mechanism if we make the kinematical assumption that the electric fleld E 
is just a function of r we obtain a flat rapidity distribution. We can also prove that the 
distribution in fluid rapidity is the same as the distribution in particle rapidity. We will also 
obtain renormalized expressions for e(r) and p(t) (non-equilibrium dynamical equation of 
state) . 

In the pair production problem we have shown that the interpolating phase space 
number density is given by the Bogoliubov function (7.7) : 

dN/drjdkr,dk±dx± = |/3(A;^, kj_, r) p (9.10) 

we are interested in transforming from dr]dkrj to dz dy where y is the particle rapidity 
y = l/2\n[{E^ + k,^)/{E^ - k,^)]. 

One can show that the transformation from {rj, r) to (z,t) is a canonical one (in the 
sense of Poisson brackets {ry, /c^} = {r, 0} = 1) with canonical momentum 

fc^ = —Ez + tp = —Tm± sinh(77 — y) 

Q = {Et — pz)/t) = m±cosh{rj — y) (9-11) 
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The phase space is unchanged by this change of variable thus 



dPN/{dr]dkr,dk±dx±) = d^N/dzdkgdk±dx± 

= JdN^/dzdydk_Ldx_L 



(9.12) 



where J ^ = dkri/dydrj/dz. At fixed r one can show that |J| = dz/dkn which leads to 
desired result, assumed by Landau that along the breakup surface r = constant: 



Schwinger's pair production mechanism leads to an Energy Momentum tensor which 
is diagonal in the(T, ry, x±) coordinate system which is thus a comoving one. In that system 
one has: 



Wc see in a 3 dimensional problem, the field theory in this approximation has two 
separate pressures, one in the longitudinal direction and one in the transverse direction 
and thus differs from the thermal equilibruim case. However, for a one-dimensional fiow 
we have that the energy in a bin of fiuid rapidity is just: 



which is just the (1 + 1) dimensional hydrodynamical result of (9.9). This result does not 
depend on any assumptions of thermalization. 

In the field theory calculation the expectation value of the stress tensor must be renor- 
malized since the electric field undergoes charge renormalization. We can also determine 
the two pressures and the energy density as a function of r. Explicitly we have in the 
fermion case. 



dN/dy = dN/df]. 



(9.13) 



T''" = diagonal {€{t),pi^{t),p^{t),p^{t)} 



(9.14) 




(9.15) 




(9.16) 



where 



RrAk) 

P\\{t)t^ 




g^\f-\^) -u^ {pI + m2)(7r + eAf/{Su\^) 



(9.17) 
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where 



R^^{k) = 2|/+|2 - {2u)-\u + A,7r)-i 
+ 57rA^(7r + ei)V(16c^V) 

and 

P±{r) =< Tyy >=< T,, >= (4r)-iE, J [dk]{pl{pl + m'y^Rrr - 2A7rpii?^^} + EI/2. 

(9.18) 

Thus we are able to numericaUy determine the dynamical equation of state pi = pi{e) 
as a function of r. 

10. DISCUSSION OF NUMERICAL RESULTS 

The physical quantities that we determine numerically are the time evolution of E(t), 
A(t) , and j(t). We will display in the figures the plasma oscillations and the time scale 
for field energy to be essentially transferred into pair production. The other quantities 
of physical interest we determine are the spectra of produced particles dN/dk and the 
dynamical equation of state. For comparison we have also solved the phenomenological 
transport theory with and without the quantum correction due to statistics (i.e. Pauli 
Blocking and Bose Enhancement) . In making plots for the spatially homogeneous case we 
use the dimensionless variables 0: i? = [eE/m?) A = eA/m;mt = r When the Electric 
field i? is > 1 then it is quite easy for pairs to be produced and in that regime the final 
result is independent of the initial data. We can see the approach to the tunneling regime 
by comparing in the regime. 5 < i?o < 2 the behavior of E(r). This is shown in fig. 3 for 
Eq = .5, 1, 2. Once Eq > 2 the behavior of E{t) is only weakly dependent on Eq. Once 
the pairs are produced one sees that there are plasma oscillations superimposed on which 
the electric field degrades. These figures are from early simulations for scalar QED in 1 + 
1 dimensions [Q. 

In fig. 4 we show A{t),e{t), < j{t) > for Eq = 2 for scalar QED in 1+1 dimensions 

i- 

In fig. 5 we show E and j for scalar QED in 1+1 dimensions for Eq = 4. We compare 
the naive Vlasov approach (dashed line) and the improved Vlasov approach ( dot- dashed 
line) . We notice that including Bose enhancement corrections is quite important. We also 
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notice that j max = 2 e pc so that particles continue to be produced when E is near a 
maximum. In fig. 6 We show the exact particle spectrum as well as the momentum space 
smoothed result which is compared to the Vlasov Equation. Here Eq =1 and we have 
scalar QED in 1+1 dimensions. 

In fig. 7 we show the results for E and j for i?o = 4 in QED in 1+1 dimension compared 
to the uncorrected Vlasov equation. We notice the dismal agreement. In fig. 8 we see 
the same curves compared to the improved Vlasov equation. In fig. 9 we show the exact 
spectrum of produced pairs for QED in 1+1 dimensions for Eq = 4. We notice that n(k) <1 
which expresses the Pauli Principle. In fig. 10 we compare the binned version of the field 
theory result with both the Naive and Improved transport theory. Next we present recent 
results for Scalar QED in 1+1 dimensions using boost invariant Kinematics. In fig. 11 we 
plot E,A and j vs. u=log(r) for Eo('Uo) = 4 in the boost invariant case where E is a function 
of the proper time r (not to be confused with the previous r). In fig. 12 we compare E(u) 
and j(u) with the boost invariant transport theory with and without Bose enhancement. 
Finally we present preliminary results pO[] for scalar QED in 3+1 dimensions. In fig. 13 
we show the time evolution for E(t) and j(t) in 3+1 dimensions and compare with the 
Boltzmann- Vlasov model with and without Bose-enhancement. 
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